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Abstract 

Two well-known methods for the design of quasicrystal models are used to create novel nonlinear 
optical devices. These devices are useful for efficient three- wave mixing of several different processes, 
and therefore offer greater flexibility with respect to the more common periodic nonlinear photonic 
crystals. We demonstrate applications for polarization switching as well as multi-wavelength and 
multi-directional frequency doubling. The generalized dual grid method is proven to be efficient 
for designing photonic quasicrystals for one-dimcnsional collinear devices as well as elaborate two- 
dimensional multi-directional devices. The cut and project method is physically realized by sending 
finite-width optical beams at an irrational angle through a periodic two-dimensional nonlinear 
photonic crystal. This enables two simultaneous collinear optical processes that can be varied by 
changing the angle of the beams. 

1 Introduction 

Twenty five years have passed since the 1982 discovery of quasicrystals [1], and we have yet to hnd 
a satisfying application that takes advantage of their unique combination of physical properties Q]- 
Nevertheless, interesting applications are starting to emerge that take advantage of quasiperiodic long- 
range order in metamaterials, or artificially constructed quasicrystalfl Most applications are based on 
linear photonic crystals, where quasiperiodic modulations of the index of refraction of a material are 
used in order to engineer its optical response. In particular, the fact that there are no restrictions on 
the order of the rotational symmetry of a quasicrystal is used to obtain nearly- isotropic photonic band 
gaps 0,0]. Here we focus on metamaterials in the nonlinear optical domain, where recent technological 
progress has enabled to modulate the second-order nonlinear susceptibility with micron-scale resolution 
in various materials, such as ferroelectrics (our focus here), semiconductors, and polymers. In these 
nonlinear photonic crystals the modulation can be achieved by planar techniques, thereby offering 
either one or two dimensions for modulation. Moreover, there are no photonic bandgaps in these 
metamaterials, because the first-order susceptibility, and hence the refractive index, remain constant. 
The advantage of using quasicrystals in this case is not in their arbitrarily-high symmetry, but rather 



'We refer the reader to Ref. 0| for a precise definition of the term 'quasicrystal' 
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in the fact that there is no restriction on the combinations of wave vectors that may appear in their 
reciprocal lattices (provided that the symmetry of the quasicrystal is not of particular importance 



The novel optical devices described below are based on materials that facilitate the nonlinear 
interaction between light waves in the form of three-wave mixing. These are processes in which two 
incoming waves of frequencies uJi and uj2 interact through the quadratic dielectric tensor x^"^^ of the 
material to produce a third wave of frequency 0^3 = a;i it W2; or the opposite processes in which a 
single wave spontaneously breaks up into two. Three-wave mixing is severely constrained in dispersive 
materials, where ^(k) is not a linear function, because the interacting photons must also conserve 
their total momentum. Even the slightest wave-vector mismatch Ak = ki it k2 — ks appears as an 
oscillating phase that averages out the outgoing wave, giving rise to the so-called "phase-matching 
problem." We have recently explained how one could fully solve the most general phase-matching 
problem using well-known ideas from the theory of quasicrystals [sl. The solution is based on the 
idea that in crystalj^l, whether periodic or not, continuous translation symmetry is broken. As a 
consequence, momentum conservation is replaced by the less-restrictive conservation law of crystal 
momentum whereby momentum need only be conserved to within a wave vector from the reciprocal 
lattice of the crystal. The fabrication of an efficient frequency-conversion device is therefore a matter 
of reciprocal-lattice engineering — designing an artificial crystal, from the quadratic dielectric field of 
the material x^'^\^)-, whose reciprocal lattice contains any desired set of mismatch wave vectors Ak^-'^ 
j = \ . . . N ^ required for phase matching any arbitrary combination of N three-wave mixing processes. 
In fact, the field amplitude of the output beam, in each frequency-conversion process, is linearly 
proportional to the amplitudes of each of the input beams, as well as the Fourier coefficient of the 
relevant mismatch wave vector (lo| . 

The idea of using a one-dimensional periodic modulation of the relevant component of the quadratic 
dielectric tensor, for the purpose of phase matching a single three- wave process, was suggested already 
in the early 1960's [ill , [l^ . 13 1, and is termed "quasi-phase matching". Since then this appro ach has 



approacn nas 



been generalized using more elaborate one-dimensional , [3] and two-dimensional 

designs, but only as ad hoc solutions for multiple processes. We argue that engineering the reciprocal 
lattice, of a nonlinear photonic quasicrystal, to contain any desired set of mismatch vectors — a task that 
25 years of research in quasicrystals have taught us how to solve — provides the most general solution 
for the long-standing problem of multiple phase-matching. Here we describe a number of novel optical 
devices that have actually been fabricated using these ideas, and tested experimentally. In Sec. [5] we 
describe devices that have been designed using the dual-grid method in order to engineer the required 
nonlinear photonic quasicrystals. These devices attest to the general nature of the quasicrystal-based 
solution to the multiple phase-matching problem. In Sec. [3] we show that for collinear devices, in 
which all the participating waves propagate in the same direction, a certain degree of flexibility can be 
obtained by using a physical realization of the cut-and-project approach, in which one generates the 
required 1-dimensional photonic quasicrystal by cutting through a periodic crystal in 2 dimensions. 



^We refer the reader to Ref. [§| for a detailed discussion on "What is a crystal?" 
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2 Reciprocal-lattice engineering with the dual-grid method 



After selecting the nonlinear medium of choice and the appropriate operating temperature for a 
given frequency-conversion application, one can calculate the required set of mismatch vectors Ak^-'), 
j = 1 . . . N , that are to be engineered into the reciprocal lattice of the nonlinear photonic quasicrystal. 
As we have described in Ref. 0], once these mismatch vectors are known, one can use de-Bruijn's 



dual grid method [20], as generalized by Gahler and Rhyner [211] and later by Rabson, Ho, and 
Mermin [22, 23], in order to design the required quasicrystal. If the N mismatch vectors are integrally- 
independent, one simply uses each vector Ak^-'^ j = 1...N, to define a family of equally-spaced 
parallel lines, separated by 27r/]Ak(-5'^ ], and oriented in the direction of Ak^-'). The set of all N 
families constitutes the dual grid, which is then used in the standard manner j^] to define a set of N 
tiling vectors a^-') , and to calculate the integral linear combinations of these tiling vectors that form the 
vertices of the tiles in the desired quasicrystal. The final step is to decorate each tile with an optimal 
motif, i.e. to decide which regions of the tile will be altered such that the relevant components of 
the quadratic dielectric tensor x^^-* are positive, leaving the remaining background with an unchanged 
negative x^^^- If N mismatch vectors happen to be linearly-dependent, one has to consider the 
pros and cons of selecting a linearly-independent subset vs. using the full linearly-dependent set for 
generating the dual grid [a]. 

2.1 One dimensional implementation 

As a rather simple demonstration of our general solution to the phase-matching problem, we have 



recently implemented a three- wave doubler 25|]. The reader who is interested in implementing his or 
her own devices, based on our solution to phase-matching problem, is kindly referred to this article 
for a detailed pedagogical explanation of each and every step. This is a one-dimensional device that 
is able to phase match three collinear second harmonic generation processes simultaneously, taking 
three input beams with wavelengths 1530nm, 1550nm, and 1570nm, and producing three output 
beams at twice their frequencies. We note that the three processes implemented in this device are 
independent, and therefore could also be phase matched by fabricating a sequence of three periodic 
nonlinear photonic crystals, one for each frequency doubling process. Nevertheless, we were able to 
show that it is more efficient to simultaneously phase match all three process in a single quasiperiodic 
nonlinear photonic crystal. 

Here we implement a more elaborate one-dimensional device utilizing our method, which is a 



nonlinear polarization switch 26|]. This device takes a beam of frequency u and polarization y, 



propagating in the x direction, and converts some of its energy into a beam of the same frequency. 



propagating in the same direction, but with polarization z. It is a cascaded device |27l ] in which the 
output of the first process is used as the input of the second process. The first is a standard (type I) 
second harmonic generation process with ujy + ujy ^ 2ujy, followed by a (type II) difference frequency 
generation process, in which the lower-frequency waves have different polarizations, 2ujy — uJy iOz- 

We implement the polarization switch using a LiNbOs ferroelectric crystal at a temperature of 
150° centigrade. Because linear dispersion is different for the differently-polarized beams, the mis- 
match vectors for the two processes are not equal, and are calculated using tabulated properties of 
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Figure 1: (Color online) The polarization switch device, (a) The Fourier transform of the device, 
showing strong Bragg peaks at the desired mismatch vectors, 0.39/im~^ and 0.66^m~^. (b) A section 
of the one-dimensional nonlinear photonic quasicrystal. The two colors represent the two different 
values, positive and negative, of the nonlinear polarization. The smallest element size is 5.35/um, 
corresponding to one half of the 10.7/im-wide strips. 
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LiNbOs [23, |29(] to be 0.39/im-i for the first process, and 0.66^m ^ for the second process. Using a 
one-dimensional version of the dual-grid method 25|] we calculate the sequence of two tiling vectors, 
of lengths 10.70/xm and 7.63/zm, that generates a one-dimensional quasicrystal with the desired wave 
vectors in its reciprocal lattice. The actual one-dimensional quasicrystal is realized by using strips 
of widths 10.70/im and 7.63^m, arranged in the x direction according to the quasiperiodic sequence 
determined by the dual-grid method. We use numerical optimization to find the optimal duty cycles 
for the two strips. We find that the best efficiencies for the desired processes are obtained when using 
duty cycles of 0.5 and 1 for the 10.7/im and 7.63/im tiling vectors respectively. This means that one 
half of each 10.7//m strip is positively-poled and the other half negatively-poled, while the narrow 
7.63//m strips are completely positively-poled. A section of the resulting one-dimensional photonic 
quasicrystal is shown in Fig. [IJb) . For a single process the best efficiency is achieved with a periodic 
nonlinear photonic crystal with a 0.5 duty cycle for its single tile. In that case the relevant Fourier 
coefficient has the value of 2/7r ~ 0.636. This value is comparable to twice the Fourier coefficients at 
the required mismatch vectors of the polarization-switch device, shown in Fig. Wi^)- However here we 
phase-match two processes simultaneously. 



2.2 Two dimensional implementation 

Next we demonstrate the versatility of our approach by designing a multi-directional second har- 
monic generator. This is a two-dimensional non-collinear device that takes three beams of wavelength 
1550nm, propagating at angles of 0° and ±20°, and generates three output beams at the same direc- 
tions at twice the input frequency. Using Stoichiometric LiTaOs and operating at 100° centigrade, 
the magnitude of the mismatch vectors, for all three processes, is calculated [30(] to be 0.301^m~^. 
The required two-dimensional photonic quasicrystal must clearly be symmetric with respect to the 
propagation direction of the input beam. Using a two-dimensional version of the dual-grid method 0], 
we find that it is generated by three tiling vectors, given in polar coordinates by a^^^ = (7.54/im, Z0°), 
a(2) = (31.31//m, Z76.9°), and a^^) = (31.31/im, Z - 76.9°). This yields two types of tiles— a rhombus 
whose edges are 31.31/im long, and a parallelogram whose edges are 31.31//m and 7.54^m long, the 
latter appearing in two mirror-related orientations. 

For ease of fabrication we limit ourselves to decorating the center of each tile with a positively- 
polarized circle, leaving the remaining background negatively-polarized. We employ numerical opti- 
mization only for determining the radii of the circles for the different tiles. We find that the best 
efficiencies for the desired processes, with the largest magnitudes of the relevant Fourier coefficients, 
are obtained for using a maximally-inscribed circle within the rhombic tile, and using no circles within 
the parallelograms. An image of the decorated tiling, as designed by our procedure, is shown in 
Fig. ^a). The calculated spectrum of the device is shown in Fig. [2||b). The reciprocal lattice vectors 
that phase match the three processes are indicated in this figure. The calculated magnitudes of the 
Fourier coefficients for the three mismatch vectors are 0.1 and 0.19 for the 0° and ±20° processes 
respectively. 
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Figure 2: (Color online) The multi-directional second harmonic generator. (a)Image of the device made 
of circular motifs, superimposed with the underlying quasiperiodic tiling, (b) Calculated diffraction 
diagram. The arrows indicate the reciprocal lattice vectors used to phase match the three second 
harmonic generation processes. 
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3 Physical realization of the cut and project method 



We have recently described an alternative scheme applicable to collinear devices 3l|]j in which one 



generates the required 1-dimensional photonic quasicrystal by a process that could be thought of as 



a physical realization of the cut-and-project method [3^, l33|]. Here we wish to elucidate some of its 



geometric features, and by doing so to emphasize its advantage for producing tunable devices. The 



reader, interested in actual implementation details of this scheme, is kindly referred to Ref. 311] • The 
basic idea is to fabricate a two-dimensional periodic crystal and employ the cut-and-project method to 
obtain a one-dimensional quasiperiodic crystal, capable of phase matching two independent collinear 
frequency-conversion processes. The cut is realized by taking advantage of the fact that the input 
beam is not an idealized plane wave of infinite transverse extent but actually has a finite width W, 
for example with a Gaussian profile. Thus, the interaction of the beam with the nonlinear medium is 
restricted to a strip-like region of width W along the propagation direction of the beam. Only those 
lattice sites of the two-dimensional crystal that fall within this strip in the transverse, or perpendicular 
(_L), direction contribute to the phase matching, and are effectively "projected" onto a one-dimensional 
quasicrystal along the propagation, or parallel (||), direction. 

The Fourier transform of the original 2-dimensional periodic crystal contains Bragg peaks at wave 
vectors K that form a periodic reciprocal lattice. As we know from the cut-and-project method, the 
Fourier transform of the strip-like interaction region is 1-dimensional. Each 2-dimensional Bragg peak 
at wave vector K gives rise to a 1-dimensional Bragg peak at the parallel component K|| of the original 
wave vector, whose intensity depends on the perpendicular component of the same 2-dimensional 
wave vector. Because the Fourier transform of a Gaussian is also a Gaussian, if for example the beam 
has a Gaussian profile and we ignore any spreading of this profile as it propagates, then the dependence 
on K_|_ is a simple Gaussian. Thus, given a pair of collinear mismatch wave vectors, one simply needs 
to find the appropriate angle with which to cut through the two-dimensional structure, so as to obtain 
two parallel projections K|| with the required mismatch values, preferably of wave vectors K with a 
small Kx component. Thus, even with a prefabricated two-dimensional periodic crystal, one has the 
ability to tune the device by varying the cut angle through the crystal, enabling the use of one device 
for different combinations of frequency conversion processes. 

To illustrate how this approach can be used, we consider a nonlinear photonic crystal built upon 
an oblique periodic lattice, defined by the primitive vectors a^^^ = (6.2/xm, Z — 75°) and a^^^ = 
(7.4/im, Z17°), by associating a positively-poled circular motif of radius 2.5/im with every lattice 
point. We wish to phase match two collinear second harmonic generation processes of fundamental 
beams with wavelengths 1550nm and 1047. 5nm. If we use stoichiometric LiTaOa and operate at 
room temperature the phase mismatch values are calculated to be 30|] Aki = 0.297/zm~^ and Ak2 = 



0.820//m~^, respectively. By varying the operating temperature we can change the required mismatch 
wave vectors, and by selecting the propagation direction of the input beams we can vary the projected 
K|| components. A simulation of all possible combinations of temperature and propagation angles 
that satisfy the phase matching conditions for the two processes is shown in Fig. [3l The simulation 
was carried out for an interaction length of 1mm and a 20//m-wide square-shaped beam profile. Darker 
shades correspond to higher efficiencies. Each parabola corresponds to a given reciprocal lattice vector 
K of the two-dimensional crystal. The apex of the parabola corresponds to phase-matching using the 
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Figure 3: Phase matching using the cut-and-project method. The calculated combinations of temper- 
ature and propagation direction that satisfy the phase-matching requirements for two collinear second 
harmonic generation processes using the oblique crystal, described in the text. The working point 
conditions are denoted by the intersection of the two straight lines at input angle a = 0° and tem- 
perature T = 23°C. The two simulated cases are: (a) for an input beam of wave length A = 1550nm, 
where phase matching is realized by a projection of the (1,0) reciprocal lattice vector; and (b) for an 
input beam of wave length A = 1047.5nm, where phase matching is realized by a projection of the 
(0,1) reciprocal lattice vector. 
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whole reciprocal lattice vector (K^ = 0), while the other parabola points correspond to projection- 
based phase matching (K^ / 0). The intersection of the straight lines is the desired working point of 
this device at room temperature, where both processes are phase-matched simultaneously. 
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